Ablowitz-Ladik system with discrete potential. I. 
Extended resolvent 



A. K. Pogrebkov*and M. C. Prati f 
February 7, 2008 



Abstract 

Ablowitz-Ladik linear system with range of potential equal to {0, 1} is con- 
sidered. The extended resolvent operator of this system is constructed and the 
singularities of this operator are analyzed in detail. 



1 Introduction 

Our aim in this article is to study the spectral theory of the matrix operator L(w), 

L m,n( W ) = S m,n-1 ~ ( J [j w \ S m>n , (1.1) 

m,n G Z, w G C, 

every element of which is a 2 x 2 matrix, S m n is the Kronecker symbol and we omitted 
a 2 x 2 unit matrix factor in the term S m<n -x. Our attention is concentrated to the case 
where values of both potentials, r„ and s n , are equal to and 1: 

r„,s„G{0,l}, neZ. (1.2) 

Moreover, we consider here the case of potentials with finite support, i.e., for every given 
potential there exist finite k and K, k < K, k, K G Z — lower and upper borders of the 
support — such that 

r n = s n = 0, n < k - 1, n>K+l. (1.3) 
The corresponding linear problem, 

L(w)$ = 0, (1.4) 

is the Ablowitz-Ladik problem [1,2] which is known to be a discretized version of the 
Zakharov-Shabat linear problem. And like the latter the Ablowitz-Ladik problem is 
associated to a variety of differential-difference integrable equations, such as discrete 



mKdV equation, difference KdV, Toda chain, etc., [3]. Problem (1.4) describes also 
discrete systems with nonanalytic dispersion relations [4] . 
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The Ablowitz-Ladik problem is also known [5,6] to be associated to difference- 
difference nonlinear equations, that are related to some class of cellular automata, i.e., 
dynamical systems in a discrete space-time with values belonging to some finite field, 
say, Fa- Cellular automata attract great interest in the literature because of the wide 
range of their applications in different sciences, from physics to biology, from chemistry 
to social sciences. Detailed references for these applications can be found in [7-11]. These 
automata are also subject to intensive mathematical study, sec for example [12-22]. It 
is just thi s kin d of applications of problem ( |l.4[ ) that motivated our specific choice of 



condition (1.2) on potential. 

The problem of the investigation of (|l.4[ ) by means of the inverse scattering transform, 
as it was performed in [3], becomes obvious if we write down this equation explicitly: 

*Ma=(£ n = 0,±l,±2,.... (1.5) 

In the standard approach to the study of the spectral problems, the main objects of 
the theory — the Jost solutions — are determined by their asymptotics at n — ► +oo an d 



n — > — oo. A solution given by its asymptotics at n — ► — oo is swept from the left by (1.5) 



But in order to construct the Jost solution given by its asymptotics at n — > +oo, one 



has to invert the matrix in the r.h.s. of (L5). The determinant of this matrix is equal to 



1 — r n s n , so in the standard appr oach the condition r n s n ^ 1 must be fulfilled. In the 



case where the potential satisfies (1.2) this means that for every n either r n or s n must 



be equal to zero [6]. Such condition drastically restricts the class of potentials of the 



type (1.2), so our aim in this and forthcoming publications is to elaborate an extension 



of the inverse scattering transform method to the case where both r„ and s n can be 



equal to 1. Let us also emphasize that, imposing condition (1.2) on the potentials, we 
do not use here the condition r n ,s n £ ¥2- As was speculated in [23] the problem of the 
integrability of the cellular automata or, more precisely, the problem of existence of the 
Lax representations must be solved in terms of the exact equalities, and not in terms of 
equalities on some finite field. 

The fact that some matrix(-matrix) operator L is analogous to a differential one 
is reflected in the property that matrix elements L m ^ n are different from zero only for 



uniformly bounded values of |m — n\. In the case of (1.1) we have 1 > m — n > 0. 
Consequently, we can apply the resolvent approach [24], [25] to investigation of the 
Ablowitz-Ladik problem. The preliminary results of our investigation were published 
in [26]. 

The resolvent approach is based on the following extension of the operator L{w): 

L m Jw,h) = h n ~ m L m Jw), (1.6) 



where h is a real non- negative parameter. In particular for the operator (1.1) we have 

L m ,n(™> h ) = ^m,n-l ~ u n{w)5 m ^ (1.7) 

where we introduced 

i/:)^ + L° r s)' neZ ' (L8) 



and a is the Pauli matrix 173, 



1 

-1 
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If we have some infinite matrix-matrix operator A m ^ n (h) depending on a parameter 
h we can associate to it the Laurent series 



+00 



A((,(',h)= ]T C m C n A m Jh), CC'eC, |CI = IC'l = i. (1-9) 



m,n— — 00 



In what follows we consider matrices A m n (h) such that the series (1.9) are convergent in 
the sense of Schwartz distributions in £, C (|C| = IC'I = 1) an( l h (h>0). The elements 
A m ,n(h) are reconstructed by means of the formula 

A m , n {h)= f j ^^ A(CC,h). (1.10) 

ICI=i IC'l=i 

In order to explain the meaning of the extension ( |l.6| ) let us introduce the function 
(distribution) 

A((,z) = A(Ce largz ,e larg Mz|) (1.11) 



where C, z S C, |C| = 1; by (\L§ and flLJ) 

A(c,z) = ^z n " m r m 4 re , n - (i-i2) 

m,n 

Then the above mentioned similarity of matrix and differential operators means that 
^4(£, z) depends on z and z^ 1 polynomially. Let us mention that if we have two objects 
of this kind, A and B, their product (composition) is defined as follows: 

+00 

(AB) m>n (h) = Yl A m,iih)B l!n {h), (1.13) 
/— — 00 

(AB)(t,C',h)= j 2m A((,C,h)B(C,C',h), (1.14) 
IC l=i 

(AB)(C,z)= j ^L A ((C,zC)B(C,z), (1.15) 
IC'l=i 



where the left hand sides of these equations are related through of (1.9)— ( 1. 12). The 
main object of our investigation is inverse M(w) of the operator L(w) extended by (|L 



L(w)M(w) = I, M(w)L(w) = I. (1.16) 



In matrix notations the first equality thanks to (1.7) has the form 

hM m+lin (w, h) = 8 m n + u m (w)M m n (w, h). (1.17) 
In order to define this inversion in a unique way we introduce 

Definition 1. A solution M(w) of ( 1.16| ) is called extended resolvent of the operator 
L{w) if M(w, h) is a Schwartz distribution with respect to £ and £' and a sectionally 
continuous function of h, h > 0. 

Let us first consider the case of zero potential, i.e., r n = s n = 0. Then the resolvent 
which we denote by Mq(w) obeys the following equation 

hM +1 (w, h) = 5 + w a M ^ n (w, h). (1.18) 
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It is convenient to rewrite this equation using representation (1.12) 

(Cz~w")M (w 1 C,z)^S c (C-l) 1 
where we introduced the S- function on |£| = 1, 

oo 

*c(c-i)= E c n > 



n— — oo 



so that 



^o(C-i)/(C) = /(i) 



I c I = i 



for an arbitrary test function /(C) on the contour. Then 

M (w,(;z) = (Cz-w^SdC-l), 



so that by ( |l.ll| ) 



MoK £ C, h) = (Cft - tO _1 * c (C/C - i), 



or by QL.1CD 

M , m ,„(«;, /i) = | - \)6{m > n + 1) 

-0(K|-/i)0(n > m)}, 
where we introduced the matrices 

0(/*-KI) = 



0(\w a \-h) 



0(h-\w\) 

o <?(/i-i/H) 

6(\w\-h) 

o e(i/H-/i) 



(1.19) 

(1.20) 
(1.21) 

(1.22) 
(1.23) 

(1.24) 



(1.25) 

Here we have to make some comments. First, by ( |l.6[ ), all expressions h m ~ n L m . n (w, h) 
are independent on h and equal to L m>n {w), see ([Ll|). On the contrary h m ~ n M . m , n (w, h) 
essentially depends on h and it is just this dependence that guaranties that Mo(w, £ £') h) 
exists as a distribution in £ £. Second, any solution of the homogeneous equation 
Lq(w)Mq(w) = is pr oportional to S(h — |ui CT |), where the matrix <5-function is defined 
in analogy with ( 1.25 ). Thus we see that the condition set on Mq(w, £ £, h) in Defi- 
nition 1 to be a sectionally continuous function of h enables us to define the resolvent 
Mq (w) uniquely. In what follows we consider the case of a nontrivial potential satisfying 
conditions (jl.2|) and (|l.3|). 



2 Extended resolvent of the regularized operator 



The specific problem connected with equation (1.17) is, as was mentioned above in the 
discussion of Eq. ( |l.4| ), that if r n = s n = 1 the matrix u n (w) is not invertible. Thus, first 
of all we have to introduce some regularization of u n (w), say, 



u n {w) -» u n (w,X) 



(Ar„s„ + l)w r n 
s n 1/w 



= A(l - det u n (w))— +u n (w). 



(2.1) 
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This substitution regularizes singular only u n (i.e., such that detu„ = 0), leaving all 
other u n untouched. Indeed, by (1.2) detw„ equals either ot 1. Then 



aetu n [w,A) | A detu n («>) = 0. 
Thus we start with the regularized operator 

L m ,n(w, A, h) = ft<S m ,n-l ~ u m(w, A)<5 m 



(2.2) 



(cf. (pq)), i.e., by 

L(w, A) = L(w) - AL>, 
where we introduced the diagonal operator 

1 + cr 



A, 



-(1 - det u n )5„ 



(2.3) 



(2.4) 



Correspondingly, we denote the extended resolvent of the regularized operator as M{w, A). 
It obeys equations (cf. ( 1 . 1 6|) ) 



L(w, \)M(w, A) = I, M(w, X)L(w, A) = I, 



that by means of (2.3) can be written in the form 

[L(w) - \D]M (w, A) = /, 
M(w,X)[L(w) - XD] = I. 

Properties of M(w, A) in the limit A — > are studied in the next section. 
Let for simplicity write 

M mtn = h m - n M m Jw,X,h), 



(2.5) 



(2.6) 
(2.7) 



(2.8) 



i.e. we omit for a while dependencies on w, A, and h. Then Eq. (2.6) takes the form 

M ,n+l,n = $m,n + »m W m,». ( 2 - 9 ) 

where dependence of u m on w and A is also omitted. It is easy to check that for any 
m > ml we have from ([2.9|) 



M m<n = 0{m>n+ l)B(n > m') [] u t + \ J\ uA M m , >n , 



where we introduced the notation 

6»(m > n) = 
and the ordered product of matrices, 



=n+l 



1, m>n, 
0, n > m + 1, 



(2.10) 



(2.11) 



m — 1 
Z=n+1 



U m~l U m~2 ■ ■ ■ U n+lJ m > 71, + 2, 

1, m = n + 1, 



(2.12) 
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Because of Eqs. (O), O), and (O) 



u = w 



n < k — 1 , 



n > K + 1, 



(2.13) 



i.e., u is a diagonal matrix in depen dent on the regulaxization parameter A . Let us 
consider first m < k. Then by ( 2.13 ) we can rewrite ( 2.10 ) in the form 



u~ m M„ 



(m > n + l)u 



(m > n + l)u 



We see that both sides of this equality are independent either on m, or on m'; we denote 
them as F n and thus we get 



M m „ =0(m>n + l)vJ- 



+ u m F n , 



< k. 



(2.14) 



Now we chose in (2.10) ml = k and substitute M k _ n in the r.h.s. using (2.14), then 



(m > n + l)8(n > k) TT u t 

l=n+l 



+d(k>n + l) J] 



. k — n— 1 



l=k 



l=k 



where m > k. Thus the second term also obeys condition m > n+1, so that taking (1.3) 
into account we can write 



m— 1 

(m > n + 1) J^J it 

/=n+l 



Y[u t ]u k F n , m>k. 



(2.15) 



Eqs. (2.14) and (2.15) give the general solution of (1.17) for any F n . In order to fix 
it, we use the two conditions formulated above. First of all it is necessary to guara ntee 
convergency for any n of the series I] m C~ m ^m : n = y~L (h()~ m M m , n where Eq. (Q 
was used. Let us consider first the sum from — oo to k. Using ( 2.14 ) we see tha t the 
sum of the first terms is finite due to ^-function. The sum of the second terms in ( [2.14 ) 



is equal (up to a constant factor) to 5Z m __ 00 u m (hC,)~ m F n . Thanks to (2.13) this sum 
converges iff the first (second) row of matrix F n is equal to zero when \w\ < h (l/\w\ < h, 
correspondingly) . Thus the condition of convergency of this series can be written as 



9(h - \w\°)F n = 0, 



n 6 Z, 



(2.16) 



where the matrix 9- function is defined in ( 1.25 ). Let us consider now the condition 
of co nverg ency of the series C, ~ m M m _ n a t plus infinity. For this purpose we write 
Eq. (|2 15b for m > K + 1 (see (jlj) and (|t|)) as 



K 



(n > K)u~ n + 6(K>n + 1) U - K J] 



l=n+l 



K 



Y[u t | u k F n } -6{n> m)u 



m—n—1 



(2.17) 
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By (2.8) the series J2m=K+i(h() m Mm,n must be convergent. The sum of the last 
terms is finite due to the ^-function and the sum of the first terms is convergent iff 



i\-h)l 9{n > K)u~ n + 9(K > n + l)i 



-K 



K 



,-K 



rh) ukF n) 

J=k J ) 



(2.18) 



The conditions (|1|) and fl2.18| ) deter mine F n uniquely. In order to get its explicit form, 
we have to consider the four regions of continuity of the matrices (1.25): 



h > \w\, h > 1/H, 
\w\ > h > 1/H> 
l/\w\ > h > \w\, 
\w\ > h, 1/| w; > h. 

Then M m n is constructed explicitly using Eqs. (2.14) and ( 2. 15] ). 
Let us introduce the (infinite) matrix column 



(2.19) 
(2.20) 
(2.21) 
(2.22) 



x m ( w , A) = 9{m > k + 1) J| Ul {w, A) w ka + 6{k > m)w r ' 



l=k 



and row 



A" 



y n {w, A) = 9{n > K)vj- na + 6{K > n+ l)ur Ka f[ u ; (w,A). 

Z=n+1 

In what follows we also use 

X m (w, A, h) = h- m x m {w, A), Y n (w, A, h) = h n y n {w, A). 
Let a denote a constant (i.e., independent on m and n) matrix 

K 



i(w, A) = w- K,7 Y[u l (w, \)w ka 



l=k 



Then, combining (2.14) and (2.15), we get 



M m>n = 6{m > n + 1) ]J Ul + x m F n 

l=n+l 



and instead of ( 2.18| ) we can write 

9{\w°\-h)aF n = -9{\w°\-h)y n . 



(2.23) 



(2.24) 



(2.25) 



(2.26) 



(2.27) 



(2.28) 



Let us co nsider first the region fl2.19| ). Then 9{h - \w a \) = 1 and by ( |2.16| ) F n = 0, 
so that (2.28) is satisfied identically. This means that in this region 



M mj „ = 9(m > n + 1) [] 



(2.29) 



l=n+l 
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In the region (|2.20|) 



6{h-\w°\)= 1 —^, 1 —^-F n = Q, 



where Eq. ( |2.16| ) was used and by (|2.2§|) 



F = 1 + a 



2a 



1.1 



Analogously in the region ( 2.21 ) we have that 

1 - a 



F„ 



2a. 



-Vn 



(2.30) 



(2.31) 



2.2 



Finally, in the region ( 2.22j ) 9{h — \w a \) = 0, thus Eq. (2.16) is satisfied identically 
and ( [2.18 ) takes the form 



All this enables us to write that 

< 

m— 1 

M m,n = 6(m>n+l) Y[ ui~ x m Ty n , 

l=n+l 

where we introduced the matrix T independent on m and n, 

0, h > \w\ and h > l/\w\, 

1 + cr 



(2.32) 



(2.33) 



r = < 



1 - a 
2a 2 2 



\w\ > h and h > l/\w\, 
h> \w\ and l/\w\ > h, 



(2.34) 



a , \w\ > h and l/|iw| > h. 
Let us mention that thanks to ( ^.13 ) we can rewrite Eqs. (2.24) and (2.32) as 



F n = -6(n > k) w - ka ^Il u iJ -0(k>n + l) w - {n+1)a , (2.35) 
and then after some simple calculations we get instead of ( 2.33p 

M mjn = -B(n>m)lf[u l \ (2.36) 

\l=rn J 

in the region ( 2.22| ). 

The results of the above construction can be summarized as the following 
Theorem 1. The extended resolvent M(w, A) of the i-operator (1.7) regularized 

by (2.1) exists, is unique and equals to 



M mn (w,X,h) =h n - m 9(m>n + l) ]J «iKA) 

l=n+l 

-X m (w,X)T(w,X,h)Y n (w,X), 



(2.37) 
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where X m , Y n , and T are given in Q2.23| ), ( |2.24| ), ( |2.25[ ), and fl2.34| ). 

To prove the theorem we need to check first of all that the double series 



+00 

E 



C m C n M m , n (w,X,h), 



(2.38) 



where (2.8) was taken into account, converge in the sense of Definition 1. Then it is nec- 
essary to prove th at (f2.37 ) obeys both (regularized) equations (1.16), i.e., equations ( |2.5|) . 

The solution ( 2.37| ) by construction is the unique solution of the first equation in fl2.5| ) 
for which the series 2m=-oo C~ m Mm,n{w, A, h) converge. Convergency of the series in 

as well as the s econd equation, are p roved analo gous ly. Both these equations easily 
follows from ( ^.37|) , if we notice that by fl2.23| ) and ( |2.24D 



x m +i( w , x ) = u m (w,X)x m (w,X), 



In other words x m and y n are solutions of the equation (1.5) regularized by 
dual. By means of ( 2.25 ) we can also write these equations as 



L(w, X)X(w, A) = 0, Y(w, X)L(w, A) = 0. 



(2.39) 
(2.40) 

l|) and its 



(2.41) 



We see that formally X and Y are right and left annulators of operator L. The existence 
of these annulators does not contradict ( |2.5| ), i.e., the existence of the inversion of L as 
both series J2 m (~~ m X m (w, A, h) and ^ n ( n Y n (w, A, h) are divergent, so X m and Y n do 
not belong to the space mentioned in discussion of Eq. (1.9) and in Definition 1. The use 
of such quantities can be avoided if, say, in the region \w\ > h, l/\w\ > h we use instead 
of (2.37) the equality 



M m . n = -h n - 



(n > m) I J^U;(w, A) 



(2.42) 



l—m 



that follows from (|J) and ( |2.36| ) 



3 Extended resolvent of the original operator 



In order to get the resolvent of the extended original operator ( |l.7[ ) we need to consider 
the behavior of (2.37) in the limit A — > 0. The existence of this limit depends on the 
regions ( ^19|)-(pT22^ Indeed, the only origin of a singularity in ( |2.37|) is matrix T, as 
follows from ( 2.23| ) and (2.24). Its limits in the first three regions of (|3§ exist by fl2~2l| . 
Let /<H _ 



a(w) 



lim a(w, X) 
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-Ka 



..ka 



(3.1) 



This expression is finite and nonzero for generic w. Zeroes of ai : i(w) and 02,2 (w) if they 
exist in the corresponding regions give bound states of operator L(w) and will be studied 
in the following publi catio n. In the region ( [2.22 ) a^ 1 (w, A) has pole at A = 0, as follows 
from the last line of ( 2.34 ) . 
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To describe the multiplicity of this pole we introduce q(m,n), m < n, number of the 
degenerated matrices M;(u>) on the interval [m, n], i.e., 



q(m,n) = ^(1 - det U;(u;)), 



(3.2) 



which is independent on w as detu^w) equals either 0, or 1. Let also 

Q = q(k,K). 



Then by (2.1) and (2.2€) we have that 



(3.3) 



a (w, A) = 



, — ko 



K 

n 

,l=k 



\w(l — det it;) — h Ui(w) 



,.Ko 



(3.4) 



where we introduce the matrices 

. , , 1/w —r n 



W 



r n 
a- 



n £ Z, 



(3.5) 



which are the inverses of u n (w) in the case where detu„ = 1 (cf. (!.£)). From (3.4) it 
follows that a -1 has pole of order Q at A = and we can write 



tCj')(«,) 



A-' 



+ 0(A), A^O. 



The residues are equal to 



K 



(Q-j) 



t&(w) = W Q - ka JJfijH 



, Z=fc 



(3.6) 



(3.7) 



where by definition 



0) 



n 



m<d <•••<(;,-<" \i=l / 



(3.8) 



for any m < n. 

By ( p. 23 ), (2.25), and (2.24) X m (w, A) and F„(w, A) are polynomials in A, so that we 
have the Laurent expansion 



M(ti;, A, ft) = M( w , A, fc) + J] — 



(3.9) 



3=1 



\w\ > h and l/|w| > h, 
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where M is the regular part of the series. The residues can be calculated explicitly 
by Eqs. (2.37) and (3.6), but in order to work with objects belonging to the space 
mentioned in Definition 1 it is reasonable to use the representation (2.42). Then by 
means of notations (3.2) and (|3.8|) we get 



(3.10) 



M$ n (w, h) = -h n - m d(n > m)9(q(m,n) > j) x 

(q{m,n)-j) 



XW 



q(m,n)—j 



Now the resolvent M(w, h) of the original (extended) L-operator (1.7) can be defined 



as 



M(w, h) 



M(w,X,h) | A=0 , 
M(w,X,h) | A=0 , 



h > \w\ or h > l/\w\, 
\w\ > h and l/\w\ > h. 



(3.11) 



Let us consider region \w\ > h, l/\w\ > h in detail. Inserting the expression (3.9) 
into Eqs. ( |2.6D and (2.7) and using (3.11) we derive that 

L{w)M{w) = I + DM w (w), M{w)L{w) = I + M ( - 1) (w)D, (3.12) 

L(w)M u) (w) = I + DM u+1 \w), (3.13) 

M<- j) (w)L(w) =1 + M {]+1) {w)D, j = l,...,Q-l, (3.14) 

L(w)M ( ' Q) (w) = 0, M {Q) (w)L(w) = 0. (3.15) 

Thus we see that in this region equations ( |3.12| ) defining the resolvent are modified in 
comparison with the standard Eqs. (1.16). In [26] it was shown that a solution of (1.16) 
does not exist in this region. In order to study the properties of the residues M™ we 
can use Hilbert identity 

M(w, X')[L(w, A') - L(w, \)}M(w, A) = M(w, A) - M(w, A') 



that follows from (2.5). Taking into account (2.6) and (2.7) we can rewrite it in the form 



M(w, A) - M(w, A') = (A' - X)M(w, X')DM(w, A). 



(3.16) 



(3.17) 
(3.18) 



Substituting M(w, A') as in Eq. (3.9), we get in the limit A' — > that 

M(w, A) - M(w) = [XM(w) - M {1) (w)]DM(w, A), 
M (j \w) = [M u+1 \w) - XM^(w)}DM(w,X), 

where j = 1, . . . , Q, and we put by definition 

M<> Q+1 \w) = 0. (3.19) 

Now we insert the expansion ( 3.11 ) in the Eqs. ( |3.17 ) and ( |3.18 ), and by ( 3.17 ) in the 
limit A — > we derive that 



M^\w)DM(w) = M(w)DM {1) (w), 

AfW(w) = M{w)DM {3+1 \w) - mW(w)DMW(w), 

j = l,...,Q. 



Then from (3.18) we have 

MW(«j) = M u+1 ' ) (w)DM(w) - MW(w)DMW(w) 



(3.20) 
(3.21) 



(3.22) 
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(symmetric to (3.21)) and 

M^ +1 \w)DM^{w) = M^{w)DM^ l+1 \w), 

M (j) (w)DM(w,X) = \m u+1 \w)D[M(w,X) - M(w)], 
A 

j,l = 1,...,Q, 



(3.23) 
(3.24) 



where ( |3.1l|) and ( |3.19[ ) were used. By fl3.22| ) we have MW'(id) = -M { - Q \w)DM i - 1 \w) 
and by $.2j ) M^\w)DM(w,X) = 0, so that also M^\w)DM{w) = by (^ll|) . 
Then for j = Q-l in ( ^22]) we have that M^^\w) = -M^- 1 \w)DM ( - 1 \w) and 
from (|3.24|) that 

AfW-^^DM^, A) = 0, so that again 1 \w)DM(w) = 0. Continuing in this way 

we prove that 



and 



M {l) {w) = -M {l) {w)DM {1) {w), 
M {l \w) = ~M {1 \w)DM {3) {w), 



M W (w)DM(w) = 0, M{w)DM {l) (w) = 0, 
j = 1, . . . ,Q, 



(3.25) 
(3.26) 



(3.27) 



where (3.26) and the second set of equations in ( [3. 27 ) is derived in an analogous way 
from (3.16), if A and A' are interchanged. 
Now by ( 3.23 ) we get 



M ij \w)DM il \w) = M^- k \w)DM il+k \w), 



(3.28) 



where j, I = 1, . . . , Q, k > 0, j — k > 1, I + k < Q + 1. If j and I are such that j — k > 1 
and / + k = Q + 1, i.e., I + j > Q + 2, then the r.h.s. of ( [3.28]) is equal to zero thanks 
to ( |3.19 ) . On the other side if / + j < Q + 1 we can chose in ( |3.28| ) fc = j — 1 and then 
by ( ]3.25| ) or ( |3~26| ) we get finally 



M {3 \w)DM {l) (w) = -9{Q + l>l+ j)M ( '+^ 1} (u;), 
j,l = l,...,Q, 



(3.29) 



that replaces relations (3.23), (3.25), and (3.26). 

This concludes the construction of the resolvent M m ^ n (w,h) of the extended L- 
operator (1.7). As we have seen, this resolvent is discontinuous at |iu| = h and \w\ = 1/h. 
In a forthcoming paper we show that study of this discontinuity leads us to modifica- 
tion of the Jost solutions and spectral data, corresponding to the case of the discrete 
potential (|L3). 
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